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NATIONAL FOREWORD 

This Indian Standard (Part 2) which is identical with ISO 1 1 843-2 : 2000 'Capability of detection — Part 2: 
Methodology in the linear calibration case' issued by the International Organization for Standarization (ISO) 
was adopted by the Bureau of Indian Standards on the recommendation of the Statistical Methods for 
Ouality and Reliability Sectional Committee and approval of the Management and Systems Division Council. 

An ideal requirement for the capability of detection with respect to a selected state variable would be that the 
actual state of every observed system can be classified with certainty as either equal to or different from its 
basic state. However, due to systematic and random distortions, this ideal requirement cannot be satisfied 
because: 

— in reality all reference states, including the basic state, are never known in terms of the state 
variable. Hence, all states can only be correctly characterized in terms of differences from basic 
state, that is, in terms of the net state variable. 

— in practice, reference states are very often assumed to be known with respect to the state variable. 
In other words, the value of the state variable for the basic state is set to zero; for instance in 
analytical chemistry, the unknown concentration or the amount of analyte in the blank material 
usually is assumed to be zero and values of the net concentration or amount are reported in terms 
of supposed concentrations or amounts. In chemical trace analysis especially it is only possible to 
estimate concentration or amount differences with respect to available blank material. In order to 
prevent erroneous decisions, it is generally recommended to report differences from the basic state 
only, that is, data in terms of the net state variable. 

— the calibrations and the processes of sampling and sample preparation add random variation to the 
measurement results. 

In this standard, the following two requirements were chosen: 

— the probability is a of detecting (erroneously) that a system is not in the basic state when it is in the 
basic state; 

— the probability is |3 of (erroneously) not detecting that a system, for which the value of the net state 
variable is equal to the minimum detectable value (x^), is not in the basic state. 

The text of ISO Standard has been approved as suitable for publication as an Indian Standard without 
deviations. Certain conventions are, however, not identical to those used in Indian Standards. Attention is 
particularly drawn to the following: 

a) Wherever the words 'International Standard' appear referring to this standard, they should be read as 
'Indian Standard'. 

b) Comma (,) has been used as a decimal marker in the International Standard while in Indian Standards, 
the current practice is to use a point (.) as the decimal marker. 
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Indian Standard 
CAPABILITY OF DETECTION 

PART 2 METHODOLOGY IN THE LINEAR CALIBRATION CASE 



1 Scope 

This part of ISO 1 1 843 specifies basic metlnods to: 

— design experiments for tine estimation of tine critical value of the net state variable, the critical value of the 
response variable and the minimum detectable value of the net state variable, 

— estimate these characteristics from experimental data for the cases in which the calibration function is linear 
and the standard deviation is either constant or linearly related to the net state variable. 

The methods described in this part of ISO 11 843 are applicable to various situations such as checking the 
existence of a certain substance in a material, the emission of energy from samples or plants, or the geometric 
change in static systems under distortion. 

Critical values can be derived from an actual measurement series so as to assess the unknown states of systems 
included in the series, whereas the minimum detectable value of the net state variable as a characteristic of the 
measurement method serves for the selection of appropriate measurement processes. In order to characterize a 
measurement process, a laboratory or the measurement method, the minimum detectable value can be stated if 
appropriate data are available for each relevant level, i.e. a measurement series, a measurement process, a 
laboratory or a measurement method. The minimum detectable values may be different for a measurement series, 
a measurement process, a laboratory or the measurement method. 

ISO 11843 is applicable to quantities measured on scales that are fundamentally continuous. It is applicable to 
measurement processes and types of measurement equipment where the functional relationship between the 
expected value of the response variable and the value of the state variable is described by a calibration function. If 
the response variable or the state variable is a vectorial quantity the methods of ISO 11 843 are applicable 
separately to the components of the vectors or functions of the components. 

2 Normative references 

The following normative documents contain provisions which, through reference in this text, constitute provisions of 
this part of ISO 1 1843. For dated references, subsequent amendments to, or revisions of, any of these publications 
do not apply. However, parties to agreements based on this part of ISO 11843 are encouraged to investigate the 
possibility of applying the most recent editions of the normative documents indicated below. For undated 
references, the latest edition of the normative document referred to applies. Members of ISO and lEC maintain 
registers of currently valid International Standards. 

ISO 3534-1 :1 993, Statistics — Vocabulary and symbols — Part 1: Probability and general statistical terms. 

ISO 3534-2:1 993, Statistics — Vocabulary and symbols — Part 2: Statistical quality control. 

ISO 3534-3:1 999, Statistics — Vocabulary and symbols — Part 3: Design of experiments. 
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ISO 1 1 095:1 996, Linear calibration using reference materials. 

ISO 1 1843-1 :1 997, Capability of detection — Part 1: Terms and definitions. 

ISO Guide 30:1 992, Terms and definitions used in connection with reference materials. 

3 Terms and definitions 

For the purposes of this part of ISO 11 843, the terms and definitions of ISO 3534 (all parts), ISO Guide 30, 
ISO 11095 and ISO 11843-1 apply. 

4 Experimental design 

4.1 General 

The procedure for determining values of an unknown actual state includes sampling, preparation and the 
measurement itself. As every step of this procedure may produce distortion, it is essential to apply the same 
procedure for characterizing, for use in the preparation and determination of the values of the unknown actual 
state, for all reference states and for the basic state used for calibration. 

For the purpose of determining differences between the values characterizing one or more unknown actual states 
and the basic state, it is necessary to choose an experimental design suited for comparison. The experimental units 
of such an experiment are obtained from the actual states to be measured and all reference states used for 
calibration. An ideal design would keep constant all factors known to influence the outcome and control of unknown 
factors by providing a randomized order to prepare and perform the measurements. 

In reality it may be difficult to proceed in such a way, as the preparations and determination of the values of the 
states involved are performed consecutively over a period of time. However, in order to detect major biases 
changing with time, it is strongly recommended to perform one half of the calibration before and one half after the 
measurement of the unknown states. However, this is only possible if the size of the measurement series is known 
in advance and if there is sufficient time to follow this approach. If it is not possible to control all influencing factors, 
conditional statements containing all unproven assumptions shall be presented. 

Many measurement methods require a chemical or physical treatment of the sample prior to the measurement 
itself. Both of these steps of the measurement procedure add variation to the measurement results. If it is required 
to repeat measurements the repetition consists in a full repetition of the preparation and the measurement. 
However, in many situations the measurement procedure is not repeated fully, in particular not all of the 
preparational steps are repeated for each measurement; see note in 5.2.1. 

4.2 Choice of reference states 

The range of values of the net state variable spanned by the reference states should include 

— the value zero of the net state variable, i.e. in analytical chemistry a sample of the blank material, and 

— at least one value close to that suggested by a priori information on the minimum detectable value; if this 
requirement is not fulfilled, the calibration experiment should be repeated with other values of the net state 
variable, as appropriate. 

The reference states should be chosen so that the values of the net state variable (including log-scaled values) are 
approximately equidistant in the range between the smallest and largest value. 

In cases in which the reference states are represented by preparations of reference materials their composition 
should be as close as possible to the composition of the material to be measured. 
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4.3 Choice of the number of reference states, /, and the (numbers of) replications of procedure, 

J, K and L 

The choice of reference states, number of preparations and replicate measurements shall be as follows: 

— the number of reference states / used in the calibration experiment shall be at least 3; however, 7=5 is 
recommended; 

— the number of preparations for each reference state J (including the basic state) should be identical; at least 
two preparations {J = 2) are recommended; 

— the number of preparations for the actual state K should be identical to the number J of preparations for each 
reference state; 

— the number of repeated measurements performed per preparation L shall be identical; at least two repeated 
measurements (L = 2) are recommended. 

NOTE The formulae for the critical values and the minimum detectable value in clause 5 are only valid under the 
assumption that the number of repeated measurements per preparation Is identical for all measurements of reference states 
and actual states. 

As the variations and cost due to the preparation usually will be much higher than those due to the measurement, 
the optimal choice of J, Kan6 L may be derived from an optimization of constraints regarding variation and costs. 

5 The critical values y^ and x^ and the minimum detectable value x^ of a measurement 
series 

5.1 Basic assumptions 

The following procedures for the computation of the critical values and the minimum detectable value are based on 
the assumptions of ISO 1 1 095. The methods of ISO 1 1 095 are used with one generalization; see 5.3. 

Basic assumptions of ISO 1 1095 are that 

— the calibration function is linear, 

— measurements of the response variable of all preparations and reference states are assumed to be 
independent and normally distributed with standard deviation referred to as "residual standard deviation", 

— the residual standard deviation is either a constant, i.e. it does not depend on the values of the net state 
variable [case 1], or it forms a linear function of the values of the net state variable [case 2]. 

The decision regarding the applicability of this part of ISO 11 843 and the choice of one of these two cases should 
be based on prior knowledge and a visual examination of the data. 
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5.2 Case 1 — Constant standard deviation 

5.2.1 Model 

The following model is based on assumptions of linearity of the calibration function and of constant standard 
deviation and is given by: 

Yfj a bxi eij (1) 

where 

X; is the symbol for the net state variable in state /; 

s^j are random variables which describe the random component of sampling, preparation and measurement 
error. 

It is assumed that the ij are independent and normally distributed with expectation zero and the theoretical 
residual standard deviation cr:f,y ~ A^(0;o-^ . Therefore, values Yij of the response variable are random variables 
with the expectation £'(3^j) a bxf and the variance v(y;-^) cr^, not depending on x,-. 

NOTE In the cases in which J samples are prepared for measurement and each of them is measured L times so that J L 

measurements are performed altogether for reference state ;, then y,y refers to the average of the L measurements obtained on 
the prepared sample. 

5.2.2 Estimation of the calibration function and the residual standard deviation 

In accordance with ISO 1 1095, estimates (see note) for a, b and 0-2 are given by: 
/ J 



{Xi-x){yij-y) 
My 1 



(2) 



a y bx (3) 

/ J 



a^ ' 



I J 2 . . . . 

' 1; 1 



[jij a bxi (4) 



The symbols used here and elsewhere in this part of ISO 1 1 843 are defined in annex A. 

NOTE Estimates are denoted by a symbol '^ to differentiate them from the parameters themselves which are unknown. 

5.2.3 Computation of critical values 

The critical value of the response variable is given by: 



1 1 J^ 



\ K I J i „„ 
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The critical value of the net state variable is given by: 



ah 1 x^ 



Xc to^95M- - -— (6) 

b \K I J s,.,, 

^o,95( ) is the 95 %-quantile of the f-distribution with =/• 7 -2 degrees of freedom. 

The derivation of these formulae is given in annex B. 

5.2.4 Computation of the minimum detectable value 

The minimum detectable value is given by: 



.0-1 1 x^ 



b \K I J s^^ 

where 

( ; ; ) is the value of the noncentrality parameter determined in such a way that a random variable 
following the noncentral f-distribution with =7-7-2 degrees of freedom and the noncentrality parameter 
, t{v\ , satisfies the equation: 

where ty ( v) is the (1 ct:)-quantile of the r-distribution with v degrees of freedom. 
The derivation of this formula is given in annex B. 
For a = p and v 3, a good approximation for 5 is given by 

( ; ; ) 2fi ( ) (8) 

if v = 4 and a = p = 0,05, the relative error of this approximation is 5%; t^ {v) is the (1 «)-quantile of the 
f-distribution with v=lJ 2 degrees of freedom. 

Table 1 presents d{v\ a; P)\or a = p= 0,05 and various values of v. 

For a = p and v 3, x^ is approximated by 



1 1 J^ 



^d 2fo,95(v')yJ- -r-: 2^0 (9) 

b\K I J s„ 
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Table 1 — Values of the noncentrality parameter for a= ;g= 0,05 and vdegrees of freedom 
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S( \a;P) 
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8( \a;P) 


V 


5{ \a;P) 


2 


5,516 


19 


3,415 


36 


3,354 


3 


4,456 


20 


3,408 


37 


3,352 


4 


4,067 


21 


3,402 


38 


3,350 


5 


3,870 


22 


3,397 


39 


3,349 


6 


3,752 


23 


3,392 


40 


3,347 
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3,673 


24 


3,387 


41 


3,346 


8 


3,617 


25 


3,383 


42 


3,344 
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3,575 


26 


3,380 


43 


3,343 


10 


3,543 


27 


3,376 


44 


3,342 


11 


3,517 


28 


3,373 


45 


3,341 


12 


3,496 


29 


3,370 


46 


3,339 


13 


3,479 


30 


3,367 


47 


3,338 


14 


3,464 


31 


3,365 


48 


3,337 


15 


3,451 


32 


3,362 


49 


3,336 


16 


3,440 


33 


3,360 


50 


3,335 


17 


3,431 


34 


3,358 






18 


3,422 


35 


3,356 







5.3 Case 2 — Standard deviation linearly dependent on the net state variable 

5.3.1 Model 

The following model is based on the assumptions that the calibration function is linear and that the standard 
deviation is linearly dependent on the net state variable and is given by: 



Yij a bxi ij 



(10) 



where 



Xj, a, b and Y-j are as defined in 5.2.1 and the s^j are independent and normally distributed with expectation 
E[£^■| = and variance: 

V{8ij) a^ixi) (c dxi ^ (11) 

i.e., the residual standard deviation is linearly dependent on x 

(t{xi) c dxj (12) 

The parameters of the model, a, b, c and d are estimated in a two part procedure as given in 5.3.2 and 5.3.3. 



5.3.2 Estimation of the linear relationship between the residual standard deviation and the net state 
variable 

The parameters c and d are estimated by a linear regression analysis with the standard deviations: 




(13) 
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as values of the dependent variable S and with the net state variable x as the independent variable. Since the 
variance V(S) is proportional to a^, a weighted regression analysis (see references [1] and [2] of the Bibliography) 
has to be performed with the weights: 

Wi -^— ^ (14) 

a (jc,) (c cfa,) 

However, the variances o^{x^ depend on the unknown parameters c and d that have yet to be estimated. 
Therefore, the following iteration procedure with weights: 

Wqi ^ (15) 

\ qii 

is proposed. At the first iteration, {q = 0), cto, = i,-, where the .?; values are the empirical standard deviations. For 
successive iterations q = ^,2, ... 

^qi Cg dqXi (16) 

calculate with the auxiliary values: 
/ 

Tq 1,1 W^i\ 

i 1 
/ 

Tq 1,2 WqiXi ; 

1 1 

/ 

Tq 1,3 WqiXi^' (17) 

i 1 

/ 

Tq 1,4 Wqi^i'' 

i 1 

/ 

Tq 1,5 WgiXiSi 

i 1 

and 

Tq XzTq 1,4 Tq ^2Tq 1,5 , , 

(^q 1 2 ('°' 

^q XVq 1,3 ^9 1,2 

and 

r^ 1,1?;? 1,5 Tg X2Tq 1,4 

"g 1 o (ly) 

?? 1,1^9 1,3 '^q 1,2 

This procedures converges rapidly so that the result for ^ = 3; 

0-3 Cg fifgX ; 
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can be considered, with 3 (x), c^ q an6 d^ t/, as tine final result: 

^{x)=\+d{x) (20) 

5.3.3 Estimation of the calibration function 

The parameters a and b are estimated by a weighted linear regression analysis (see references [1] and [2] in the 
Bibliography) with the y,-^- as values of the dependent variable, x, as values of the independent variable and 
weights: 

1 



Hxi) 



where 

^{xj) is the predicted value of the variance at x- according to equation (20) 
with: 





I 


T^ 


J wr, 




i 1 




/ 


T2 


J WiXi ■ 




i 1 


T3 


/ 
J Wjxf; 




i 1 




/ J 


Ta 


J Wijij ; 




M ;• / 




/ J 


T5 


J WjXiyij 




i \j I 


the estimates for a and h are 


a 


T3T4 T2T5 

T{r^ Ti 


b 


T{Tz T2T4 
T{T^ Ti 



5.3.4 Computation of critical values 

The critical value of the response variable is given by: 



(21) 



(22) 



(23) 



Jc=«+fO,95(v)J-^+ ^ + ^ '' (24) 



8 
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and the critical value of the net state variable is given by: 

(25) 



Xc = 


f0,95(^) 


-2 . -2 


-2 


e 


T2IT, 







(26) 
73 TilT^ (27) 



.2 1 /_ . .-^ "^ 



/■7-2 . . . . 



W, 



\yij -a- bxi (28) 



and fo,95('^) 's the 95 %-quantile of the r-distribution with v= IJ -2 degrees of freedom; 5„w is defined in annex A. 

5.3.5 Computation of the minimum detectable value 

The minimum detectable value is given by: 



;,^4PM+ 1 + ^-2 (29) 

where 

S = S(v; a; /^) is the value of the noncentrality parameter as defined in 5.2.4. 
Since " ^(x^) depends on the value of x^ yet to be calculated, x^ has to be calculated iteratively. 

The iteration starts with "(x^)o "0 ^"d results in x^q; for the next iteration step ~(x^)^ "(-*^do) '^ computed and 
used in the formula for x^, resulting in x^-^,... In many cases even the first iteration step does not change the value 
of x^ appreciably; an acceptable value for x^ is obtained at the third iteration step. 

6 Minimum detectable value of the measurement method 

The minimum detectable value obtained from a particular calibration shows the capability of the calibrated 
measurement process for the respective measurement series to detect the value of the net state variable of an 
observed actual state to be different from zero, i.e. it is the smallest value of the net state variable which can be 
detected with a probability of 1 -/? as different from zero. This minimum detectable value differs for different 
calibrations. The minimum detectable values of different measurement series for 

a particular measurement process based on the same type of measurement process, 

a type of measurement process based on the same measurement method, or 

a measurement method 

can be interpreted as realizations of a random variable for which the parameters of the probability distribution can 
be considered characteristics of the measurement process, the type of measurement process or of the 
measurement method, respectively. 
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If, for a particular measurement process, m consecutive calibrations have been carried out in order to determine the 
minimum detectable value of the net state variable x^, the m minimum detectable values x^i, x^2' ■■■ ^dim can be 
used to determine a minimum detectable value of the measurement process under the following conditions: 

a) the measurement process is not changed; 

b) the distribution of the values x^j is unimodal and there are no outlying values x^; 

c) the experimental design (including the number of reference states, /, and the numbers of replications of 
procedure, J, K and L) was identical for each of the calibrations. 

Under these conditions the median of the values x^j, for / = 1, ..., m, is recommended as the minimum detectable 
value of the measurement process; if another summary statistic of the values x^^ is used instead of the median, the 
statistic used shall be reported. 

If any of these conditions are violated, the minimum detectable value of the measurement process is not sufficiently 
well-defined and the determination of a common value shall not be attempted. 

If the same measurement method is applied in p laboratories and for each of them a minimum detectable value of 
the measurement process within the laboratory were to be determined, then under the same conditions as for the 
determination of the minimum detectable value of the measurement process, the median of the p minimum 
detectable values of the laboratories is recommended as the minimum detectable value of the measurement 
method; if another summary statistic of the minimum detectable values of the laboratories is used instead of the 
median, the statistic used shall be reported. 

7 Reporting and use of results 

NOTE Examples of the determination of critical and minimal detectable values are given In annex C. 

7.1 Critical values 

For decisions regarding the investigation of actual states only the critical value of the net state variable or of the 
response variable is to be applied. These values derived from a calibration of the measurement process are 
decision limits to be used to assess the unknown states of systems included in this series. Looking at consecutive 
calibrations of the same measurement process, the critical values may vary from one calibration to another. 
However, since each of the critical values is a decision limit belonging to a particular measurement series, it is 
meaningless to calculate overall critical values across calibrations and logically inappropriate to use these as 
critical values. 

If a value of the net state variable or of the response variable is not greater than the critical value, it can be stated 
that no difference can be shown between the observed actual state and the basic state. However, due to the 
possibility of committing an error of the second kind, this value should not be construed as demonstrating that the 
observed system definitely is in its basic state. Therefore, reporting such a result as "zero" or as "smaller than the 
minimum detectable value" is not permissible. The value (and its uncertainty) should always be reported; if it does 
not exceed the critical value, the comment "not detected" should be added. 

7.2 Minimum detectable values 

The minimum detectable value derived from a particular calibration shows whether the capability of detection of the 
actual measurement process is sufficient for the intended purpose. If it is not, the number J, K or L may be 
modified. 

A minimum detectable value derived from a set of calibrations following the conditions mentioned in clause 6 may 
serve for the comparison, the choice or the judgement of different laboratories or methods, respectively. 
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Annex A 

(normative) 

Symbols and abbreviations 



a intercept in tine expression y a bx e 

a estimate of tine intercept a 

b slope in tine expression y a bx s 

S estimate of tine slope b 

c intercept in the expression (x) = c + dx for the residual standard deviation 

c estimate of the intercept c 

d slope in the expression {x) = c + dx for the residual standard deviation 

d estimate of the slope d 

E ( ) expectation (of the random variable given in the brackets) 

I number of reference states used in the calibration experiment 

/ = 1 , ..., / identifying variable of the reference states 

/ number of preparations for each reference state 

7 = 1, ..... / identifying variable of preparations for the reference- and basic state 

K number of preparations for the actual state 

k=^, ...,K identifying variable of preparations for the actual state 

L number of repeated measurements for each preparation 

/ = 1 ,..., L identifying variable of the repeated measurements per preparation 

M multiplying factor 

m number of consecutive calibrations 

A'^ number of preparations in the calibration experiment; if the number of preparations for 

each reference state is identical, then N=I-J, and the total number of measurements in 
the calibration experiment is N-L 

q = 0,1,2, ... number of the iteration step 

s empirical standard deviation 
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J (x- B^ sum of squared deviations of the chosen values of the net State variable for the reference 

' states (including the basic state) from the average 



2 weighted sum of squared deviations of the chosen values of the net state variable for the 
reference states (including the basic state) from the weighted average 



J Wi[Xi x„) 

i 1 



T auxiliary value for the weighted linear regression analysis 

y ( ) variance (of the random variable given in the brackets) 

w; weight at X; 

Wqi weight at x, in the q^'^ iteration step 

X net state variable, X = Z-zq 

X a particular value of the net state variable 

x-^, ..., xj chosen values of the net state variable X for the reference states including the basic state 

Xg critical value of the net state variable 

■^d minimum detectable value of the net state variable 

- \ average of the chosen values of the net state variable for the reference states (including 

/ ' the basic state) 

i 1 
y„ a 

X ^-^ — estimated value of the net state variable for a specific actual state 

_ ' / ' weighted average of the chosen values of the net state variable for the reference states 

x„ wiXj / Wj (including the basic state) 

M / r 1 

Y response variable 

yf, critical value of the response variable 

yiji 1^^ measurement of the/*^ preparation of the i^^ reference state 

y^i , ..., yf,i obtained values of the response variable for the fc^^ preparation of a specific actual state 

in the measurement series 

^ K L 

ya yu average of the observed values for a specific actual state 

I J L 



1 



yiji average of the measurement values y,- 



I J L. .. ., /-' ''^' 

1 ^ J M ^ 

1 ^ 
yij T ^'ji average of the measurement values of thej"^ preparation of the j"^ reference state 

^ I 1 
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J L 

yiji 



J ^ju. 



average of the measurement values of the /'^ reference state 



Jo average of the KL measurement values at x = 

Z state variable 

zq value of the state variable in the basic state 

a probability of erroneously rejecting the null hypothesis "the state under consideration is 

not different from the basic state with respect to the state variable" for each of the 
observed actual states in the measurement series for which this null hypothesis is true 
(probability of the error of the first kind) 

in the absence of specific recommendations the value a should be fixed at «= 0,05 

p probability of erroneously accepting the null hypothesis "the state under consideration is 

not different from the basic state with respect to the state variable" for each of the 
observed actual states in the measurement series for which the net state variable is 
equal to the minimum detectable value to be determined (probability of the error of the 
second kind) 

in the absence of specific recommendations the value /? should be fixed at /?= 0,05 

S non-centrality parameter of the non-central f-distribution 

E component of the response variable measurement representing the random component 

of sampling, preparation and measurement errors 

V degrees of freedom 

diff standard deviation of the difference between the average, y , and the estimated intercept, 

a 

estimate of the residual standard deviation 
qi Standard deviation at x,- in the q^^ iteration step 

' estimate of the residual standard deviation, x = 



13 



IS/ISO 11843-2:2000 



Annex B 

(informative) 

Derivation of formulae 



B.I Case 1 — Constant standard deviation 

Under the assumptions of 5.1 and in tine case of constant standard deviation, estimations of tine regression 
coefficients, a and b, are normally distributed with expectations 

E{a) a ; Ei^b) b 

and variances: 

va ^ -^ -■,vb - 

where 

^ is the variance of the residuals of the averages of the L repeated measurements for each preparation. 

If the response variable is measured KL times at the basic state (z zq,x 0) , the difference between the average 
^0 of the KL values and the estimated intercept a follows a normal distribution with expectation: 

E\yo "5) EyyQj E\aj a a 
and variance: 



2 2 2 

Viyo a) V{yo) V {a) _ ^ ^ 2 1 ^ ^ 

K I J s^x K I J s^ 

Since (yo fl) is normally distributed, the random variable 

u ^^^ 
diff 

follows the standardized normal distribution, and the inequality: 



diff 



"0,95 



holds with probability 0,95. Since %^ is unknown it can be estimated as: 

diff „ r T 
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where 

"^ is the estimated residual variance of the regression analysis that shall be used instead. The random 
variable 

diff 
follows the f-distribution with I J 2 degrees of freedom, and the inequality: 

Jo « 



diff 

or 



f0,95(^) 



1 1 x^ 



Jo « f0,95(l')^diff « f0,95(>^) ,, „ , , 

where 

fo,95( ) 's the 95 %-quantile of the f-distribution with vdegrees of freedom, holds with probability 0,95. 
The right hand side of this inequality is the critical value of the response variable. 



1 1 x^ 



V ^ 1 J s 



and the critical value of the net state variable is 



^c =^^^ = ^0,95(1^)^, - + y7 + — 
^ b\K IJ s^ 

Similar expressions describe these values when other quantiles of the f-distribution are appropriate. 

In order to determine the minimum detectable value x^ of the net state variable, it is necessary to examine the 
distribution of (y a)/"^|ff in the case where the true value x of the net state variable is identical to the minimum 
detectable value x^ of the net state variable, x = x^. It is required to detect this state with probability 1 , i.e: 

P fo,95( )| JC ^d 1 

diff 

or 

<^diff 
If X = xj, the expectation of y is: 

E{y) a bx^ 
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and therefore 




E{y a) 


bx^ 


whereas: 




V{y a 


2 

diff 


as for X = 0. 




p y « 

diff 


f0,95(^) ^ ^d 


-P ^~ 


d - bxfi + bx^ 




^ m 


y 


a bx^ bx^ 


p 


diff diff 




'diff/ diff 


-P —. tnarAv) 



f 0,95 (1^)1 



f0,95(^) 



p[nv\S) fo,95(l')]; 

since U (y a bx^)l ^jjff follows the standardized normal distribution and "^jiff/ ^jjff independent of U follows 

the distribution of ^x'^{v)lv, the random variable T{v, ) follows the noncentral r-distribution with v degrees of 

freedom and noncentrality parameter 5; ( ; ; ) for 0,05 or other appropriate value, if required is 

determined as the value of the noncentrality parameter of the noncentral f-distribution with k degrees of freedom 
that satisfies: 

p[n ; ) h ( )] 

From: 

S ^^ 

diff 

the expression: 



b b\K I J s^^ 

for the minimum detectable value of the net state variable follows. 
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For a prognosis, the estimates of b and a are inserted into tine formula so tinat tine minimum detectable value is 
given by: 



id S 



a- 1 1 X 

b\K I J s. 



The critical value of the response variable y^ is the sum of a and a multiple of , and the critical value of the net 
state variable is a multiple of ^ jb. If, according to the recommendations, the values of the net state variable of the 
reference states are equidistantly spaced with the smallest value zero, a= 0,05 and either 

K= 1 (one preparation for the measurement of the actual state) or; 

K = J (number of preparations for the measurement of the actual state equal to this for the reference states); 

the multiplier: 



^ fo,95(i')i 



'1 1 X' 

K I J s. 



in the expressions for the critical values is a function of the number of reference states, /, and the number of 
preparations of each reference state, J, only. For some cases M is given in Table B.1 . 







Table B.I 


— Determination of the multiplier factor, M 




?ox K=^ 


I 


J 


IJ 




f0,95( ) 


M 




3 


1 


3 


1,35 


6,31 


8,52 


3 


2 


6 


1,19 


2,13 


2,54 


5 


1 


5 


1,26 


2,35 


2,97 


5 


2 


10 


1,14 


1,86 


2,12 


5 


4 


20 


1,07 


1,73 


1,86 


?ox K = J 


I 


J 


IJ 




f0,95( ) 


M 


1/ 1 x^ 


3 


1 


3 


1,35 


6,31 


8,54 


3 


2 


6 


0,96 


2,13 


2,04 


5 


1 


5 


1,26 


2,35 


2,97 


5 


2 


10 


0,89 


1,86 


1,66 


5 


4 


20 


0,63 


1,73 


1,09 



17 



IS/ISO 11843-2:2000 



B.2 Case 2 — Standard deviation linearly dependent on the net state variable 

Under the assumptions of 5.1 and in tine case of tine standard deviation being linearly dependent on the net state 
variable, the estimations of the regression coefficients, a and b, are normally distributed with expectations: 

E a a ; E\b b 

and variances: 



V a 



V b 



Tz 2 1 ^w 



2 

xxw 



rr 2 



T-^T^ Ti 



where 

^ is defined so that w,- ^ is the variance of the residuals of the averages of the L repeated measurements 
for preparation /. 

If the response variable is measured KL times at the basic state Z zq, X , the difference between the 
average y of the KL values and the estimated intercept a follows a normal distribution with expectation: 

E y a E y E a a a 

and variance: 

2 ^ -2 

Vy a Vy Va — — ^iiff 

"^ ■'1 ^xxw 

diff is unknown, but can be estimated as follows: 

-2 -2 ^ -2 

-2 _ 0,v^_ 0, T , ^w -2 

diff-^r + ^« -^r+ :^ + 

where 

"o is taken from equation (20) and '^ is the estimated residual variance of the weighted regression 
analysis, which shall be used instead. 

In analogy to case 1 the critical value of the response variable is: 



U 2 . -2 

3'c=« + f0,95(^)"diff =« + fO,95(l')J-^+ — + ^^ "^ 



and the critical value of the net state variable is: 



. _, txA diff _fa95M , 1 --^.y ^2 

Similar expressions describe these values when other quantiles of the f-distribution are appropriate. 
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These formulae include the case of constant standard deviation for which all the weights are equal to one, w,- 1 
for / 1,...,/ so that r^ I J , x„, x, i„^^ s^ and '§ '^■ 

The minimum detectable value of the net state variable is: 

b 
where, forx = x^, 

diff.xd ^(y o\x Xd) y{y\x x^ Y{a) 

For a prognosis, the estimates of h and \\\^x^ , b and: 

"- ^ 1 ^xxw 

are inserted into the formula so that the minimum detectable value of the net state variable is given by: 



Xd= — .1 + — + - 



Since ^(x^j depends on the value of x^ yet to be calculated the iterative procedure of 5.3.5 has to be used. 
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Annex C 

(informative) 

Examples 



C.1 Example 1 

The mercury content, expressed in ng/g'') of plant materials, was measured by atomic absorption spectroscopy. 
Each sample was decomposed using a microwave (MLS-1200) technique and taken up in nitric acid / potassium 
dichromate solution. These solutions were examined through a Varian VGA-76 cold vapour reduction system 
leading to a gold-plated foil concentration system (MCA-90) prior to replicated atomic absorption measurements. In 
order to estimate the calibration function and to determine the capability of detection, each of six reference samples 
representing the blank concentration (x = 0) and the net concentrations x = 0,2 ng/g; 0,5 ng/g; 1,0 ng/g; 2,0 ng/g; 
3,0 ng/g was prepared three times and each prepared sample measured once. Hence, / = 6; 7 = 3; L = 1 . 

It was assumed that the assumptions of linearity of the calibration function, constant standard deviation and normal 
distribution of the response variable hold; a and /?had been fixed in advance at 0,05. For the determination 

of the concentration of mercury in the material to be analysed, two different approaches were taken into 
consideration: 

a) one measurement would be carried out (/: = L = 1 ) ; or 

b) three samples would be prepared for measurement and each of them measured once (K = 3; L = 1) and the 
average y^ of the observed values used as the measurement result. 

The results of the calibration experiment are given in Table C.1 . 



Table C.1 — Results of the calibration experiment for the determination of mercury content 

in food or drugs 



Reference 
sample 


Net concentration of 
mercury 




Absorbance 




i 


ng/g 




yij 




1 





0,003 


-0,001 


0,002 


2 


0,2 


0,004 


0,005 


0,005 


3 


0,5 


0,011 


0,011 


0,012 


4 


1,0 


0,023 


0,023 


0,023 


5 


2,0 


0,048 


0,047 


0,048 


6 


3,0 


0,071 


0,072 


0,072 



1) 1 part per billion (ppb) = 10 9 g/g = 1 ng/g. The use of ppb is deprecated. 
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The statistical analysis yields: 
X 1,116 7ng/g 
i„ 20,425 



a 9,995 9 10 ^ 



b 0,023 74 
' 1,109 9 10^ 
Since v=N-2 = 16; 

fo,95( ) fo,95(16) 1,746; 

( ; ; ) (16; 0,05; 0,05) 3,440; 
(2fo,95( ) 3,492 
The results for the approach a) are 

critical value of the response variable [see equation (5)] y^ 0,003 05 

critical value of the net concentration [see equation (6)] x^ 0,086 ng/g 

minimum detectable net concentration [see equation (7)] x^ 0,1 73 ng/g 

the smallest absorbance value which can be interpreted as coming from a sample with a net mercury 
concentration larger than the blank concentration is y^ 0,003 05 , the critical value of the response variable; 

the smallest net concentration of mercury in a sample which can be distinguished (with a probability of 
1 0,95) from the blank concentration is x^ 0,1 73 ng/g, the minimum detectable value of the net 

concentration. 

The results for approach b) are: 

critical value of the response variable [see equation (5)] y^ 0,002 30 

critical value of the net concentration [see equation (6)] x^ 0,055 ng / g 

minimum detectable net concentration [see equation (7)] x^ 0,1 1 ng / g 

C.2 Example 22) 

The amount of toluene in 100 |il of extracts was measured using gas chromatography interfaced with a mass 
spectrometric detector (GC/MS). 100 |il samples were injected into the GC/MS system. Six reference samples 
were used and contained toluene in known amounts in the range 4,6 pg/100 |il to 15 000 pg/100 |il. Each sample 
was injected and measured four times (/ = 6, 7 = 4, L = 1 , A' = 24). The measurement results are given in Table C.2. 



2) D.M. RocKE and S. Lorenzato. A Two-Component Model for Measurement Error in Analytical Chemistry. Technometrics, 
1995, 37, pp. 181-182. 
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A look at the graphical representation of the measurement results shows that the relationship between the toluene 
amount and the response variable (peak area) is satisfactorily linear; the standard deviation of the peak area is 
linearly dependent on the amount of toluene. Under the additional assumption of normal distribution of the 
response variable the capability of detection can be determined according to 5.3. 

Table C.2 — Results of the calibration experiment for the toluene amount in 100 I extract 



(1) 


(2) 




(3) 




(4) 


(5) 


(6) 


(7) 


Reference 
sample 


Net amount 
of toluene 




Peak area 




Empirical 
standard 
deviation 


Predicted standard deviation of 
iteration 














1 


2 


3 


/ 


pg/1 00 nl 




yij 




'i 


1/ 


"2; 


"3; 


1 


4,6 


29,80 


16,85 


16,68 


19,52 


6,20 


4,56 


5,17 


5,15 


2 


23 


44,60 


48,13 


42,27 


34,78 


5,65 


7,07 


7,93 


7,92 


3 


116 


207,70 


222,40 


172,88 


207,51 


21,02 


19,73 


21,87 


21,88 


4 


580 


894,67 


821,30 


773,40 


936,93 


73,19 


82,91 


91,43 


91,57 


5 


3 000 


5 350,65 


4 942,63 


4 315,79 


3 879,28 


652,98 


412,46 


454,22 


455,02 


6 


15 000 


20 718,14 


24 781,61 


22 405,76 


24 863,91 


2 005,02 


2 046,54 


2 253,14 


2 257,23 



In the estimation procedure for c and d an iteratively reweighted linear regression analysis according to 5.3.2 is 
carried out which produces the following estimated linear regression functions: 



iteration 1 : 



'1; =3,933 23 + 0,136 174 jc,- 



iteration 2: " 2; = 4,482 84 + 0,1 49 91 1 x^ 
iteration 3: " 3,- = 4,462 28 + 0,1 50 1 85 x,- 



The corresponding predicted standard deviations are given in columns (5) to (7) of Table C.2. After the third 
iteration the results are stable so that the equation of iteration 3 can be used as the final result of part 1 of the 
estimation procedure, i.e.: 

'(;c) = 4,462 28 +0,150185 x 

^ 4,462 28 

The parameters a and b of the calibration function are estimated by a weighted linear regression analysis according 
to clause 5.3.3 with the yjj of column (3) as values of the dependent variable, x,- of column (2) as values of the 
independent variable and weights: 



1 



W: 



[xi) (4,462 28 + 0,1 501 85 X;)' 
This regression analysis yields: 



T.=J w: = 0,223 306 



1 



= 15,566 9 
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s,^ = 606,224 

a =12,218 5 

b =1,527 27 



■^ =1,059 54 



V =N~2 = 22 

f0,95( ) =f0,95(22 =1,717 

Therefore for /*: = 1 , the following are obtained: 

critical value of the response variable [see equation (24)] y^ 20,82 

the critical value of the net toluene amount in 1 00 I of extract [see equation (25)] jcq 5,63 pg. 

The minimum detectable value is calculated iteratively: 

For a 0,05, ( ; ; )= (22; 0,05; 0,05) = 3,397 (see Table 1) and with ^{xd)^ ^q the first value for x^ 

[see equation (29)] is jCcio = 11,139; it follows ^{xd)^ 6,135 2 and Xdi 14,553; 

with ^(xd)2 6,647 9 iteration 2 leads jcci2 15,627 pg/1 00 land 

with ^(xd)3 6,809 2 we get finally xd JCdS 15,967 pg/1 00 I. 

The smallest peak area which can be interpreted as coming from a sample with a net toluene concentration larger 
than the blank concentration is y^ 20,82 , the critical value of the response variable. 

The smallest net amount of toluene in a sample of 100 I extract which can be distinguished (with a probability of 
1 0,95) from the blank concentration isx^ 15,97 pg/1 00 I, the minimum detectable value of the net toluene 

concentration. 
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